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Abstract 

At nuclear matter density, electrically neutral strongly interacting matter 
in weak equilibrium is made of neutrons, protons and electrons. At sufficiently 
high density, such matter is made of up, down and strange quarks in the color- 
flavor locked phase, with no electrons. As a function of increasing density 
(or, perhaps, increasing depth in a compact star) other phases may intervene 
between these two phases which are guaranteed to be present. The simplest 
possibility, however, is a single first order phase transition between CFL and 
nuclear matter. Such a transition, in space, could take place either through a 
mixed phase region or at a single sharp interface with electron-free CFL and 
electron-rich nuclear matter in stable contact. Here we construct a model for 
such an interface. It is characterized by a region of separated charge, similar 
to an inversion layer at a metal-insulator boundary. On the CFL side, the 
charged boundary layer is dominated by a condensate of negative kaons. We 
then consider the energetics of the mixed phase alternative. We find that the 
mixed phase will occur only if the nuclear-CFL surface tension is significantly 
smaller than dimensional analysis would indicate. 



1 Introduction 



It is becoming widely accepted that at asymptotically high densities, the ground 
state of QCD with three quark flavors is the color-flavor locked (CFL) phase 
Ijl], ^ 1^, 3. In this phase, there is complete spontaneous breaking of the color gauge 
symmetry, chiral symmetry and baryon number. This phase of QCD is, as we will 
see below, a transparent, insulating superfluid. Moreover, the effective coupling is 
weak and the ground state and low-energy properties can be determined by adapting 
the BCS methods used in the theory of superconductivity ||, H, |, |^, H, ||, |10], |Tl|, 
12t O, ll^ 0, H]- At less-than-asymptotic densities, the CFL phase continues to be 



the ground state for nonzero quark masses, and even for unequal masses, so long 
as the differences are not too large ||T5|, 0. This means that the CFL phase is the 
ground state for real QCD, in equilibrium with respect to the weak interaction, over 
a substantial range of densities. 

Could the cores of neutron stars, long speculated to contain quark matter, consist 
of this remarkable phase, whose properties are calculable from first principles? To 
discover the answer, we must establish observable consequences of the CFL-core 
hypothesis. A necessary step is to understand the transition(s) between ordinary 
nuclear matter and CFL quark matter occurring with increasing depth in the 
neutron star. 

The two simplest scenarios are (1) a single sharp interface between nuclear matter 
and CFL, and (2) a mixed phase region. We construct consistent semi-quantitative 
models for both. Which is more favorable depends on the surface tension of the 
interface. We will conclude that for the mixed phase to occur, the surface tension 
must be significantly less than the value suggested by naive dimensional analysis. 

We begin with a brief summary of the general properties of the CFL phase, and 
overviews of the nuclear-CFL single interface and the mixed phase. 

1.1 CFL generalities 

The CFL phase consists of equal numbers of u, d and s quarks, and so requires no 
electrons to make it an electrically neutral, macroscopically allowed bulk phase . 



The CFL pairing energy, associated with the formation of ud, us, and ds Cooper 
pairs, is maximized when all three flavors have equal number density. This equality 
is enforced and the electrical neutrality of the CFL phase is undisturbed even in the 
presence of a nonzero strange quark mass (up to some critical value). 

Although the CFL ground state breaks the color and electromagnetic gauge 
symmetries, there is an unbroken U{1)q gauge symmetry and a corresponding 
massless "rotated" photon given by a specific linear combination of the ordinary 
photon and one of the gluons [0, |18|. The CFL state is neutral with respect to Q- 
charge. From now on, whenever we discuss "charge" or "electric field" in the CFL 
phase we will always be referring to the electrodynamics of the rotated photon that 
couples to Q. 
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At temperatures which are small compared to the superconducting gap A (of 
order tens to 100 MeV), the transport and response properties of CFL quark 
matter are dominated by the lightest excitations. There is an exactly massless 
superfluid mode, associated with the spontaneous breaking of the exact baryon 
number symmetry. The lightest charged excitations are the pseudoscalar mesons, 
which are the pseudo-Nambu-Goldstone bosons associated with spontaneous chiral 
symmetry breaking At quark chemical potential /i, these have masses of 

order ^m,m„,rfA//i, of order ten MeV [g 0, |, |, [Tl|, 0, |13, 



13. The effective 



field theory which describes these light degrees of freedom is known and at high 
enough density all coefficients in it can be determined by controlled, weak-coupling 
calculations [g |, 0, |, |, |10|, |lT| 0, [g . 



Because there are no electrons, and the charged hadronic modes are gapped, Q- 
electromagnetic flux, in the form of photons or DC electric or magnetic fields, can 
penetrate zero-temperature CFL quark matter unimpeded, obeying free Maxwell 
equations. The zero-temperature CFL phase in bulk is therefore a transparent 
insulating superfiuid, and continues to be so up to temperatures of order tens of 
MeV, well above the temperature of any neutron star more than a few seconds old. 



1.2 Single nucIear-CFL interface 

The simplest possibility for a nuclear-CFL interface is that at a single quark chemical 
potential ii^ electrically neutral nuclear and CFL matter have equal pressure, and 
there is a single first order phase transition between the two phases. We will argue 
below that for nuclear matter, which has a large nonzero electron chemical potential 
/ie, to be in stable contact with electron-less CFL matter, a charged boundary layer 
must develop, extending over both sides of the interface. One aspect of this is 
a "QCD-scale" micro-boundary region, with a length scale presumably of order 
one fermi, where the microscopic structures of the different ground states somehow 
mesh. We will idealize this interface as infinitely thin, and we do not attempt to 
analyze its details. (In the end, our ignorance of the QCD-scale physics of the micro- 
boundary will prevent us from giving a completely definitive answer to the question 
of the stability of the single interface.) Our analysis of the interface focuses on 
the physics of the charged boundary layers extending over both sides of the micro- 
boundary. These are tens of fermi wide. They can be analyzed quantitatively based 
on effective theories for the low-energy degrees of freedom above the bulk ground 
states on either side of the interface. For the CFL phase, we have such a theory in 
hand, and we know that it becomes accurate asymptotically. On the nuclear side, 
there is unfortunately nothing of comparable rigor. We will use a Walecka model, 
which is well-documented and easy to implement. Eventually, more sophisticated 
descriptions of the nuclear matter side should be analyzed. 

Why is there interesting physics on the 'macroscopic' scale, at tens of fermi? 
Consider attempting to construct a sharp interface, where neutral CFL and nuclear 
matter meet at a micro-boundary of order 1 fm, on either side of which we find 
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the respective bulk phases. Due to the difference in chemical potentials, electrons 
will ffow from the nuclear side of the micro-boundary to the CFL side. This flow 
halts only after an electric field develops, as the residual net positive charge on the 
nuclear side attracts the electrons on the CFL side, keeping them from penetrating 
too far into the CFL matter. The outcome is a charge-separated interface, with a 
layer of positively charged, electron-depleted, nuclear matter on one side and a layer 
of CFL quark matter with electrons on the other, stabilized by the resulting electric 
field. The natural length scale for such an electron boundary layer is the electron 
Debye length, Ae = /i'^y^vrTia 

em ~ 10 fm. Note that on the CFL side, the electric 
field is that associated with U{1)q whereas on the nuclear side, it is that of ordinary 
electromagnetism. 

The electrons are far from being the whole story, however. It turns out that 
although they set the length scale for the thickness of the charged boundary layers, 
they do not dominate the charge density. Due to the electric field, protons on 
the nuclear side pile up near the interface, making this layer even more positive. 
Similarly, the electric field induces a condensate of negative kaons on the CFL 
side of the interface, where the kaons are the lightest possible negatively charged 
excitations. The net negative charge density on the CFL side of the interface is 
dominated by the kaons, while the positive charge density on the nuclear side is 
dominated by the protons. We show the electron, proton and kaon number density 
profiles near a model interface in Fig. ^ in 

Due to the separation of charges, the electrostatic potential is position- 
dependent. This means that although /Xg is constant across the interface, the 
"effective" electron chemical potential fif^ = (/ig + e0) is position-dependent. 
Its value is He deep in the nuclear matter and zero deep in the CFL matter. 
The bulk CFL phase remains electrically neutral even in the presence of a large 
nonzero /ig imposed by contact with the nuclear matter because of the presence of 
a compensating electrostatic potential. 

We will calculate the macroscopic density profiles at the minimal interface 
by using a Thomas-Fermi description for the fermions and a Landau- Ginzburg 
description for the kaons, and solving a self-consistent equation for the electrostatic 
potential (essentially, the Poisson equation). The calculation is reminiscent of one 
previously performed to analyze the electric field at the interface between vacuum 
and quark matter in the absence of pairing, where electrons spill out of the quark 
matter |]T9|. In our case, electrons spill into the quark matter from the nuclear 
matter, and the protons and kaons turn out to play a major role in addition. We 
present the calculation in after laying the necessary groundwork in and 
In we describe our models of the bulk CFL and nuclear phases in detail, and in 
§1^, we analyze the interface upon making the simplifying assumption that no kaon 
condensate occurs on the CFL side. 

Let us note several interesting features of the minimal interface. 

• The proton and kaon densities are large. 
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• One must take into account the fact that whereas describes an ordinary 
electric field on the nuclear side of the interface, it describes a Q field on the 
CFL side. This means that the electric field (— V0) must satisfy nontrivial 
boundary conditions at the interface, dual to those derived for magnetic fields 
in Ref. 0. 

• From the point of view of compact star physics, the most striking feature of the 
minimal interface is probably the simple fact that it introduces a discontinuity 
in the density-vs.-radius profile of such a star. For the particular choice of 
parameters we analyze, nuclear matter with baryon density of 2.1nQ and energy 
density 343 MeV/fm^ fioats on CFL matter with baryon and energy density 
both about twice as large, meeting at an interface whose boundary layers 
are only tens of fermi thick. Here, no = 0.16 fm~^ is the nuclear saturation 
density. The consequences of this density discontinuity for the properties of 
a static compact star and for the dynamics of binary inspiral warrant much 
further investigation. 

1.3 Nuclear-CFL Mixed phase region 

As Glendenning realized the bulk energetics does not favor a single interface. 
Instead, it suggests the existence of a mixed phase region, with domains of positively 
charged nuclear matter interweaving among domains of negatively charged CFL 
matter. This phenomenon can be understood from Fig. 0, a schematic phase 
diagram of QCD in the n-nl^ plane. If one neglects electromagnetism, and thus 
allows charged bulk phases, nuclear matter is stable in the lighter (yellow) region 
of Fig. |1] while CFL matter is stable in the darker (violet) region. They meet 
along a coexistence line, where the two phases have the same chemical potentials 
and pressure, but different electric charge densities. The CFL phase is electrically 
neutral on the heavy (red) line given by = 0. The nuclear phase is electrically 
neutral along the heavy (red) line through AB. 

Two possible paths from nuclear to CFL matter as a function of increasing fi are 
depicted. In the absence of electromagnetism and surface tension, the favored option 
is evidently a mixed phase made of negatively charged CFL matter and positively 
charged nuclear matter along the segment of the coexistence line from A to D. On 
this segment, positively charged nuclear matter coexists with negatively charged 
CFL matter, so for pressures in the range Pa to Pd an overall neutral mixed phase 
can be created by choosing an appropriate volume fraction of CFL relative to nuclear 
matter. We construct the mixed phase in §^ If, on the other hand. Coulomb and 
surface energies are large, then the mixed phase is disfavored. The system remains 
on the nuclear neutrality line up to B, where there is a single interface between 
nuclear matter at B and CFL matter at C. This minimal interface, which we 
construct in with its attendant charged boundary layers, occurs between phases 
with the same /Xg, = /^s = /ic? and pressure P^,. The effective chemical potential 
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Figure 1: A schematic phase diagram showing the nuclear and CFL phases in the 
plane of quark chemical potential fi and effective electron chemical potential fif^. 
Isobars are shown as thin solid lines, and neutrality lines for nuclear and CFL matter 
are thick (red) lines. Each phase is negatively charged above its neutrality line, and 
positively charged below it. Continuation onto the unfavored sheet is shown by 
broken lines. 

fif^ changes across the interface, though, as a result of the presence of the electric 
field. 

So, is the minimal interface as previously described stable? To decide between 
the single interface and mixed phase scenarios, we must consider the cost of making 
multiple interfaces, that is the surface tension. 

The surface tension of the single interface turns out to be dominated by the cost 
of the macroscopic boundary layer, rather than the "fundamental" surface tension 
<^QCD of the micro-boundary. We calculate the boundary layer surface tension in 
§1^, and find that it is quite large: about 420 MeV/fm^. This result alone does not 
preclude the existence of a mixed phase, however. Rather, it constrains the size 
of the domains in such a mixed phase to be small compared to the 10 fm length 
scale characteristic of the boundary layers that dress a single interface. If the mixed 
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phase has small enough Wigner-Seitz cells, density profiles do not vary much within 
a single cell, and the surface energy reverts to ctqcd- In §11 we show that while a 
mixed phase region is guaranteed to occur for small enough cqcd, the single interface 
constructed in §^ is more stable than the mixed phase if ctqcd > 40 MeV/fml On 
the other hand, naive dimensional analysis suggests o"qcd ~ 300 MeV/fm^. 



1.4 Other possibilities 

Throughout the remainder of the paper, we will adopt the minimal hypothesis that 
nuclear matter composed of neutrons, protons and electrons, known to exist at 
densities of order uq, and CFL matter, known to exist at asymptotically large 
densities, exhaust the content of a neutron star. We conclude this introduction, 
however, by mentioning some of the many non-minimal possibilities. 

If the transition from nuclear to quark matter occurs at too large a density, then 
interesting complications may arise on the nuclear side; if it occurs at too small a 
density, interesting complications may arise on the quark side. 

If the transition to quark matter occurs at a high enough density, it will be 
preceded by the onset of kaon condensation and/or by the onset of nonzero 
hyperon density. Either kaon condensation or the presence of hyperons will tend 
to reduce fig- If the hyperon number densities get sufficiently large that hyperon- 
nucleon pairing occurs, it may even be possible for the transition to the CFL phase to 
occur continuously p^, |15|! (Note that the transition from ordinary nuclear matter 



to the CFL phase cannot be continuous, since it involves a change in symmetry.) 

If the transition to quark matter occurs at a low density or equivalently, since the 
relevant parameter is m^/(4/xA), if the effective strange quark mass at the transition 
is too large or the gap is too small, then the transition may occur to a less symmetric 
form of quark matter, in which CFL pairing does not occur. Ordinary BCS pairing 
may still occur between quarks of two out of three flavors. There will certainly be 
electrons present in order to maintain neutrality. The quarks which cannot undergo 
BCS pairing may well undergo LOFF pairing |^, making the matter a crystalline 
color superconductor p4| , p6| . In this non-minimal scenario, it is only as the 



density increases further that a subsequent phase transition to quark matter in the 
CFL phase occurs. 



2 Bulk Neutral CFL vs. Bulk Neutral Nuclear Matter 

In this section, we construct the free energies of electrically neutral nuclear and 
CFL matter, choose parameters, and find the chemical potential at which the two 
phases have equal pressure. These provide the basis for analyzing the density profile 
of model neutron stars containing a CFL quark core and a single interface. We 
leave that analysis, and an analysis of how the location of the interface varies as 
parameters are changed, to future work. Here, for concreteness, we settle on one 
reasonable set of parameters. 
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2.1 The CFL phase. . . 

We describe the CFL phase using the free energy |T7| 



Q r'^ 3 f / \ 

f^CFL = — P'^ip- f^)dp + — P^ [ Vp'^ + ^^s - Mdp 



3AV' 

^2 



+ B, (2.1) 



where the quark number densities are Uu = = rig = {u^ + 2A^/i)/7r^ and the 
common Fermi momentum is 



2/x-^//z2 + ^^/z-g. (2.2) 

Let us discuss each term in turn. The first two terms give the free energy of the u 
and d quarks, assumed massless, and s quarks with their mass m^, in the absence 
of interactions. However, the number densities are not what they would be in 
the absence of interactions. CFL pairing forces all fiavors to have the same Fermi 
momentum and hence the same number density, as long as is not too large |T7 . 



The next term is the leading contribution (in powers of A//i) of the CFL pairing to 
the free energy 0, and is validQ whether the interaction which causes the pairing 
is treated as some phenomenological four-quark interaction or as the exchange of a 
propagating gluon §]. (The calculation of A is quite different in these two cases, 
but the contribution to the free energy is as in ( p.l| ) as long as A/fi is small.) 
The derivation of the first three terms in ^cfl is given in Ref. which uses a 
simplified two-quark model, but the generalization to the nine quarks of the CFL 
phase is straightforward. The final term is a bag constant, which we use as a simple 
phenomenological way of parametrizing the physics of confinement. 

It would certainly be possible to include additional physical effects to what 
appears in ^cfl- For example, we are neglecting the perturbative interactions 
among the quarks ||2^. The CFL pairing, which we do include, has qualitative 
effects which play a crucial role in the following. Including the perturbative effects 
would simply have the effect of increasing ^cfl at a given fi. The fact that we leave 
this out means that to obtain a reasonable phenomenology, we must choose larger 
values of B than are typically used when the perturbative effects are included 



It is worth noting that in the CFL phase the value of corresponding to a given 
pressure P = —Q depends sensitively on P and on the bag constant B, but only 
weakly on the gap A and on m^, as long as both are small compared to {B + Py^^. 
This can be seen by rewriting (|2.1| ) as 



V = ./i7r2(i? + P) + W-A2 + ... (2.3) 



^The numerical coefficient in front of this term is quantitatively valid if the CFL condensate is 
purely in the color-3, flavor-3 channel. There must, in fact, be a small additional condensate in 
the color- and flavor-symmetric channel |[ but this makes a negligible contribution to the 
free energy. 
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to lowest order in A and m^. 

The electromagnetic properties of CFL matter will be crucial in subsequent 
sections, so we review them here. In QCD with three flavors, electromagnetism 
is described by a U{1)em symmetry which is a gauged subgroup of the flavor group 
SU{3)l X SU{3)ji. In the CFL phase, each separate SU{3) in the S'f/(3)coior x 
SU{3)l X SU{3)ji symmetry of the Lagrangian is spontaneously broken, but the 
symmetry SU {3)coior+L+R associated with simultaneous color and flavor rotations 
remains unbroken. One f/(l) subgroup of this unbroken SU{3)co1ot+l+r symmetry is 
a linear combination of f/(l)EM and a [/(I) subgroup of the original color symmetry. 
Once we are alerted to this possibility, it is not difficult to identify the appropriate 



combination of the photon and one gluon which remains unbroken |T8[. It is 
generated by 

Q = Q + r]Ts (2.4) 

where Q is the conventional electromagnetic charge generator, Tg is associated with 
one of the gluons, and in the representation of the quarks. 



Q = diag(|, — |, — |) in flavor u, d, s space, 
r/Tg = diag(— |, |, |) in color r, g, b space. 



(2.5) 



By construction, the Q-charges of all the Cooper pairs in the CFL condensate are 
zero. (For example, with these conventions red up quarks pair only with green down 
or blue strange quarks, and both these pairs have Q = in sum.) The condensate is 
Q-neutral, the U{1) symmetry generated by Q is unbroken, the associated Q-photon 
will remain massless, and within the CFL phase the Q-electric and Q-magnetic fields 
satisfy Maxwell's equations. The linear combination of the photon and the eighth 
gluon which remains massless is 

Al = ' / ' I = cos aoA, + sin aoG« (2.6) 

with 7] = 1/ v^, while the orthogonal combination 

Al = ' \ / = - sin aoA, + cos a^Gl , (2.7) 

gets a mass and the corresponding X-magnetic field experiences a Meissner effect. 
The denominators arise from keeping the gauge field kinetic terms correctly 
normalized, and we have defined the angle ao which specifies the unbroken U{1) 
via 

cosao = = ■ (2.8) 

yjlfe^ _|_ g2 

The mixing angle is the analogue of the Weinberg angle in electroweak theory, 
in which the presence of the Higgs condensate causes the and the third SUi^w 
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gauge boson to mix to form the photon, A^, and the massive Z boson. At accessible 
densities the gluons are strongly coupled (5f^/(47r) ~ 1), and of course the photons 
are weakly coupled (e^/(47r) f» 1/137), so ao rje/g is small, perhaps of order 
1/20. The rotated photon consists mostly of the usual photon, with only a small 
admixture of the Tg gluon. 

All the elementary excitations in the CFL phase (the pseudo-Goldstone bosons, 
the massive vector bosons, the gapped fermions) couple to with charges which 
are integer multiples of 

— ecoscco , (2.9) 



the Q-charge of the electron, which is less than e because the electron couples only 
to the component of J^. The only massless excitation, the superfluid mode 
corresponding to spontaneous violation of baryon number, is Q-neutral. Because all 
charged excitations have nonzero mass, and there are no electrons present, the bulk 
CFL phase at low temperatures is a transparent insulator. 

In the vicinity of the interface with nuclear matter, electrons (with charge — e), 
the negative kaons (same charge) and the Q-electric field will all play a role. 

Before attempting a comparison between I^cfl and i^nucicar, we must determine 
for what values of A, and \i the CFL phase is more stable than less symmetric 
quark matter, i.e. quark matter in the absence of CFL pairing, rendered electrically 
neutral by the presence of a nonzero electron density. In the CFL phase, even 
though /ie takes on the same value as that in the nuclear matter with which it is 
in contact, the effective electron chemical potential /ig + e0 = and the electron 



density vanishes. In unpaired quark matter, weak equilibrium imposes = A* — g/^e 
and //(i = /Xs = A* + |A*e and electrical neutrality turns out to require 

where the higher order terms are suppressed by further powers of m^//i^. We can 
therefore evaluate the difference between the free energy of neutral CFL quark 
matter and neutral unpaired quark matter. We find: 

3 3 

^^CFL - ^^unpaired = + J^"^* ' (^-H) 

to lowest nonzero order in A/// and m^///^. This means that as long as 

2 

A > ^ (2.12) 

the free energy gained from CFL pairing is greater than the free energy cost of 
maintaining equal quark number densities. This criterion for the stability of the 
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CFL phase relative to that of neutral unpaired quark matter is the analogue of 
that derived in Ref. [|17| in a simplified two-quark model, although the numerical 
coefficient in ( p.l2D differs from that in the model of Ref. . 



Actually, the decisive comparison is not that between CFL pairing and no pairing 
at all. Rather one should compare the free energy of the CFL phase with that of 
some "2SC phase" wherein standard BCS pairing occurs between either up and down 
quarks only or up and strange quarks only. The quarks which do not participate in 
the 2SC condensate may form a crystalline color superconducting phase [0, ^ |26 



in which quarks with differing Fermi momenta pair with each other, or else quarks of 
the same flavor may pair among themselves to form spin-1 condensates Either 



form of secondary pairing, single-flavor or crystalline color superconducting, make 
only negligibly small corrections to the free energy. Inclusion of the 2SC condensate 
affects the competition between CFL quark matter and less symmetrically paired 
quark matter mainly by reducing the coefficient in ( p. 12 ) from 1/4 to a somewhat 



smaller value, no less than l/vl2. Since this has little impact on the present paper, 
and the analysis includes features of independent interest, we will defer it to a 
subsequent publication. 

2.2 The nuclear matter phase. . . 

We could make many of the qualitative points we wish to make if we idealized 
the nuclear matter side of the interface as noninteracting neutrons, protons and 
electrons. However, if we want to work with estimates of the numerical values of 
the number densities near the interface which are potentially illustrative, we must 
incorporate the strong interactions among the neutrons and protons in some way. 
In this initial effort, we choose to use a Walecka-type relativistic field theoretical 



N 



model in which the nucleons interact with omega, rho and sigma mesons [Q. The 
nuclear-CFL interface could and should be studied using less simplified treatments 
of the nuclear matter side. The Lagrangian we use for the nucleon sector is given 
by 

+ \d,ad^a - \mla^ - U{a) - \v,^V>^^ 

+ \<y,V^ - \r,u ■ R^" + \mlR, ■ (2.13) 

where m*^ = mj^ — QaNO' is the nucleon effective mass, which is reduced compared 
to the free nucleon mass m^v due to the scalar field a, taken to have mo- = 600 MeV. 
The vector fields corresponding to the omega and rho mesons are given by Vf^i, = 
dfjYi, — dj^Vfj,, and R^^, = d^R^ — d^R^j^ respectively. The scalar self- interaction term 
is given by 

U{a) = \mN{gaNcrf + -A9aNo)^ , (2.14) 
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where b and c are dimensionless coupling constants. \[^Ar is the nucleon field operator 
with tm the nucleon isospin operator. The five coupling constants, g^N, QujN, dpN, 
b, and c, are chosen as in Ref. to reproduce five empirical properties of nuclear 
matter at saturation density: the saturation density itself is no = 0.16 fm"^; the 
binding energy per nucleon is 16 MeV; the nuclear compression modulus is 240 
MeV; the nucleon effective mass at saturation density is O.TSm^v; and the symmetry 
energy is 32.5 MeV. 

The model is solved in the mean-field approximation, wherein only the time 
component of the meson fields have nonzero expectation values. The symbols cr, u 
and p denote sigma, omega and rho meson expectation values that minimize the 
free energy given by [pOl 



1 / f'^Fn rkpp \ 

finuclear(/U„,/ie) = ij dk k"^ {Snik) - Pn) + J dk P {ep{k) - fXp) j 

+ I {m'y - mlu' - my) + f/(o-) - , (2.15) 



where 



en{k) = \Jk'^ + mlf'^ + g^^uj - \gpNP , (2-16) 
^p{k) = ^Jk^ + ni*/ + g^N^ + \gpNp , (2.17) 

are the neutron and proton single particle energies in the mean field approximation. 
The corresponding Fermi momenta kpn and kpp-, which minimize the free energy at 
fixed baryon and electron chemical potentials, are given by solving 



:2.i8i 



^nikpn) — Pn ; 
€p(^kpp) = Pp , 

where weak equilibrium sets pp = pn — Pe, and 

mla = g„N [risikpn) + risikFp)) - ^ , 

mluj = g^N (nikpn) + n{kFp)) , (2-19) 
= \9p {n{kFp) - n^kpn)) ■ 

The nucleon number density n and scalar density for nucleons with Fermi 
momentum kp are 

kp 1.3 
- 2 '^F 



1 

nikp) = — 5" / dkk 
Jo 



1 /""-P 
risikp) = ^ dkk 
Jo 



37r2 

2 "^Af 



1 rkp * (2-20) 



2 
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Note that in Eq. ( |2.15| ), the electron contribution has also been included. 
In bulk matter, the condition dflnuciear/dfie = 0, which enforces electric charge 
neutrality, uniquely determines fie- The magnitude and the density dependence of 
the electron chemical potential is sensitive to the value of the nuclear symmetry 
energy, parametrized in this model as the strength of the isovector interaction. 



2.3 . . . and their meeting point 

The critical quark chemical potential, /ic, above which the electrically neutral 
CFL state has lower free energy than nuclear matter is determined by requiring 
^CFhifJ'c) = f^nuciear(/in = 3/ic,/ie)- (In Fig. |I|, = = fic-) Both /Xc and the 
critical pressure Pc = — ^^cFL(Aic) depend sensitively on the high density behavior 
of the nuclear EOS and on the numerical value of the bag constant. They are less 
sensitive to the value of the gap, as Eq. ( p.3| ) suggests. For the nuclear equation 
of state described in the previous section, and for CFL-phase parameters given by: 
51/^ = 190 MeV, rus = 150 MeV and A = 100 MeV, we find that /i^ = 365 MeV. 
The corresponding electron chemical potential on the nuclear side, required to ensure 
electrically neutral bulk nuclear matter, is fie = 214 MeV. The pressure at the 
interface is 34 MeV/fm^. The baryon density on the nuclear side is n'^^^'^^^ = 2.1 rig 
and on the CFL side is ra^^^ ~ 4.3 uq. In estimating n'^^ = —dflcFh/dfi, we have 
treated A as //-independent .0 The energy density on the nuclear side is 343 MeV/fm^ 
and on the CFL side is 719 MeV/fm^. Note that with the parameters just described, 
the criterion (p.l2|) is satisfied by a factor of more than six in the quark matter at 
fie- This justifies our assumption that the quark matter is in the CFL phase. It is 
only if A were significantly smaller or if the effective strange quark mass were 
significantly larger that we would find a transition from nuclear matter to a less 
symmetric form of quark matter (with 2SC and crystalline color superconducting 
condensates) followed at a larger /i by a second transition to the CFL phase. 

As an illustration of the sensitivity to B, if we reduce B^^^ to 171 MeV, CFL 
matter with density 2.7 no becomes stable at zero pressure. To give a sense of the 
effects of CFL-pairing on these bulk properties, note that if we had not included 
the effects of CFL-pairing in f2cFL, we would have found that unpaired strange 
quark matter was stable at zero pressure only for 5^/^ = 155 MeV, instead of 
171 MeV. Looked at another way, with a fixed choice of B^^'^ (for example 190 MeV 
as we use), CFL pairing lowers the free energy of the CFL phase, thus increasing 
its pressure, and therefore reduces the values of Pc and fic at which the interface 
with nuclear matter occurs, relative to previous estimates made using unpaired 
quark matter. Although a more systematic study incorporating effects like the 
perturbative interaction among quarks would be necessary before making contact 
with various phenomenological normalizations of B ES], we expect this qualitative 



^With the parameters we have chosen, the A^fi^ term of (2.1) only contributes at the level of 
A^//i^ ~ 10% to the number density, so any reasonable //-dependence of A wiU contribute below 
this level. 
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feature to be robust. As a consequence, CFL-quark matter will extend closer to the 
surface of a compact star than previously estimated for unpaired quark matter. 

Leaving construction of neutron star models and the systematic exploration of 
the dependence of the location of the interface on parameters like B, rris and A 
for future work, we turn now to the physics occurring at the CFL-nuclear interface, 
wherever it occurs. 



3 The Minimal Interface Without Kaons 

The bulk calculation of the previous section gives us /i and /Xe at the nuclear- CFL 
interface. We now set up the calculation of the number densities near the interface. 
We choose a geometry where a sharp phase boundary exists ai z = 0, the region 
z < contains nuclear matter and the region ^ > contains CFL quark matter. 
As previously explained, the physics of concern here occurs on the 10 fm distance 
scale, so we do not attempt to resolve the internal structure of the micro-boundary, 
instead encapsulating its properties in boundary conditions at 2; = 0. The region 
near the interface is characterized by a positive charge density on the nuclear side 
and a negative Q-charge density on the CFL side. As a result, there is an electric 
field E{z) on the nuclear side and a Q-electric field Eq{z) on the CFL side, which 
we express in terms of an electrostatic potential 0, 

E{z) = V0 2; < (nuclear) , 
Eq{z) = V0 z>0 (CFL) . 

The boundary condition for perpendicular electric field at the interface is simply 

Eq{0+) = cosao^(O-) , ^^2) 
i.e. dz(f){0~^) = cos ao 920(0") , 



which is dual to the condition derived for perpendicular magnetic fields in Ref. |jT8 
The underlying physics is that the electric field entering the CFL phase from the 
nuclear phase is resolved into a Q component, which penetrates into the CFL region, 
and an orthogonal X-component which is screened out. The CFL condensate carries 
X-charge, and screens X-flux on a length scale of order 1/A, which is short compared 
the electron Debye length Ag. This justifies our treating it by a boundary condition. 
In reality, we see from ( |2.8| ) that cosao is very close to unityQ so in most of our 
calculations we set cosao = 1. In Fig. ^ we illustrate the effects of ao on the 
interface using the exaggerated case cosao = 1/2. 

The boundary conditions at infinity, deep in the nuclear and CFL phases, are 

^^^(-oo) = (nuclear) , 

d,(j){+oo) = (CFL) . ^ ■ ^ 



As the Q dielectric constant in the CFL phase is shghtly different from one , the ratio in 
21) is shghtly smaller than cosaoi but is still very close to unity. 
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These follow from the fact that the star as a whole is neutral, as are the bulk CFL 
and nuclear phases. 

The electron chemical potential /le must be constant across the interface, 
otherwise electrons would flow. However, the electron density is controlled by the 
effective electron chemical potential 

fj,f{z) = /ie + e(f){z) z<0 , 

^ IJ,e + e(f){z) z> , ^ ' ^ 

which is z-dependent, enabling the electron density to vary across the interface. 
The fact that the CFL phase is neutral in the absence of electrons means that 
ljLf^{+oo) — 0, allowing the boundary conditions to be expressed as 



(f){—oo) — (nuclear) , 

0(+oo) = -He/e (CFL) . 

To obtain the density profiles we solve the Poisson equation 



(3.5) 



dz"" 



epqiz) , (3.6) 



subject to the boundary conditions above, where pg is the electric charge density. 
On the CFL side of the interface, e is replaced by e. We solve Poisson's equation 
in the local density or Thomas-Fermi approximation, where the charge density is 
written in terms of position-dependent Fermi momenta kppiz) and kpeiz) for protons 
and electrons respectively: 

Pq(^) = ^2 ■ (3-7) 

(In this section we neglect the hght charged bosons in the CFL phase, which 
give a further contribution to the charge density which cannot be described via 
a local Fermi momentum; these effects will be included in the next section.) The 
Thomas-Fermi approximation is valid if 0, and hence kpp and kpe, vary only on 
length scales which are long compared to l/kpe and l/kpp- That is, we require 
dkF^ p/dz -C ^Fgp- We have checked that the profiles we find satisfy this condition 
at the 5% level. This confirms our expectation that the local Fermi momenta vary 
on a length scale of order Ag ~ 10 fm. In addition to justifying our use of the 
Thomas-Fermi approximation, this justifies our treating the 1 fm scale physics via 
boundary conditions at 2; = 0. 

For non-interacting electrons and protons, the local Fermi momenta can be 
simply expressed in terms of the ^^-independent chemical potentials //g and and 
the electrostatic potential (t){z) 

^ ' ^ ^ (3.8) 



klpiz) = {^p-^f{z)Y- 



m 



N 
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Figure 2: Particle number density profiles and electrostatic potential near a sharp 
interface between nuclear matter and CFL quark matter. Distance is measured in 
units of the electron Debye length Ag = /i~^A/vr/4aem ~ 10 fm. In constructing 
these boundary layers, we have neglected the possibility of a kaon condensate on 
the CFL side. The transition occurs at = /in/3 = 365 MeV and an electron 
chemical potential /Zg = 214 MeV. Solid curves show results which include effects 
due to strong interactions between baryons. Dashed curves show results where they 
are neglected. 



where /ip = /i^ — /ie = A/yUg + is the proton chemical potential. It is fixed 
by requiring electrical neutrality deep in the nuclear matter region: kpp{—oo) = 
kpe{—oo). To include the effects of strong interactions among the nucleons, instead 
of using the equation above for kpp one must solve the equations of the Walecka 
model given in §0, which yield kpp and k function of /x„ and /if (Eqs. M 

and (gl9D ). 

The Thomas-Fermi equation for is an ordinary second-order differential 
equation, which we have numerically solved using both shooting methods and 
relaxation methods. Our results for the electrostatic potential and the resulting 
particle density profiles are given in Fig. ^ The distance z from the interface is given 
in units of the electron Debye screening length Ae = ^/ttJAo^. The dashed 
curves are the density profiles assuming nuclear matter made of noninteracting 
neutrons, protons and electrons. The solid curves show the effects of including the 
strong interactions among the nucleons using the Walecka model described in §0. 
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Notice that because protons and neutrons are strongly coupled, the neutron density 
readjusts slightly near the interface to ensure constant baryon chemical potential. 

Because the electrons are taken to be massless, the electrostatic potential and 
the electron density in the CFL phase exhibit power law falloff kpe = IJ-l^ ^ 1/z. 



4 The Minimal Interface With Kaons 

Because the CFL condensate breaks chiral symmetry, the lightest charged degrees 
of freedom in this phase are the resulting Nambu-Goldstone excitations. These 
are light but not massless once nonzero quark masses are taken into account. The 
pseudoscalar bosons can be analyzed in terms of an effective field theory with the 
Lagrangian |^, |^ 

= :^Tr (SoSSoS"^ + vldiJ^d'i:^) + c (detMTr(M-^S) + h.c.) , (4.1) 

where we have retained only terms with two derivatives. The color singlet 

5 = exp iv^i3//,r transforms under SU{3)l x SU{3)r and describes the octet of 
pseudoscalar mesons defined by the 3x3 matrix B. The masses of the pseudoscalar 
mesons associated with chiral symmetry breaking can be obtained from the structure 
of the mass terms in the above Lagrangian |ll|. The pion and kaon masses are 
given by 

"^^± = -F^^sirriu + rud) , m^i = |^md(m„ + m^) . (4.2) 

Thus, the kaon is lighter than the pion, by a factor of m^/ (m^ + m^) We have 
neglected instanton effects which induce a small (gg) condensate in the presence of a 
CFL condensate. The resulting additional contribution to the meson masses [|l], ||, 
is likely small |3^, |^. 

At asymptotically high densities, the full theory is weakly coupled and the 
coefficient c appearing in the mass term, the velocity of the modes and the decay 
constant can therefore all be calculated from first principles. Up to possible 
logarithmic corrections, they are given by |^, ^, ^, |TT], |T2|, |T^, 0] 

21-81og2 2 2 1 

3o7r^ 3 ZTT^ 

At the densities of interest to us, these asymptotic expressions can certainly not be 
trusted quantitatively, although we shall use them in order to be concrete. 

As first noted in the CFL context by Schafer [^, a nonzero electron chemical 



potential may drive the condensation of negatively charged Goldstone bosons. This 
renders the matter negatively charged, and therefore cannot occur in bulk CFL 
matter. The bulk CFL phase has /i^^ = 0. We find that meson condensation does 
occur within the negatively charged boundary layer on the CFL side of the interface. 
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At first encounter, one might expect that negative kaons (which have charge — e hke 
the electrons) condense if fif^ > rriK and negative pions condense if fi'f > rrin. 
This is not the case, however. Even if fif^ > m,r, it is always more favorable to 
make (lighter) kaons than pions. Simply put, pions can always decay into kaons. 
We have confirmed by explicit calculation that if we incorporate the possibility of 
both pion and kaon condensation, only kaon condensation occurs. For simplicity of 
presentation, therefore, we retain only the kaon fields in the meson matrix B.^ We 
characterize the kaon condensate by an amplitude 6 and rewrite charged kaon fields 
as 

/ A^+\ 

B = \ with = ^ exp [T^f^ft] . (4.4) 

\K- V 2 



Substituting this into Eq. ( |4.1| ), we find that free energy density associated with the 
K~ condensate is given by 

^]^(/i, /if) = -| ((/if )2 sin^ e + vliyef + Am\ sin\\e)) . (4.5) 

Varying VLk with respect to (which enters in VLk via /if) yields a new contribution 
to the right hand side of the Poisson equation. Varying VLk with respect to 9 
yields a new differential equation because of the term. It turns out, however, 
that wherever the kaon density is appreciable, this spatial gradient term in VL^ is 
negligible. More precisely, for the 80% of kaons closest to the interface, the 
term is always 5% or less of the kaon free energy. We therefore drop the spatial 
gradient term, and solve explicitly for the condensate amplitude by minimizing VLk, 
obtaining 6{z) from cos6'(2;) = m|'/(/if (z))^. Thus, we can rewrite the free energy 
in the following compact form: 

M...f)^-i/.Vf^(l-2^ + ^). (4.6) 

The kaon number density is then 



nK{^^, /.f ) = ""T'^^ ' = fl /^f 1 - 7^ • (4-7) 



gjj^ _ i^"KVP)Pe ) _ f2 ,,eff / i "J^ 

,efi 

All these expressions are only valid for /if > mj^; for smaller /if, there is no kaon 
condensate: d = rix = 0. 

When rriK "C /ie the Debye screening length of kaons is given by Xk = 
(\/47raem /vr)^^ ^ 10 fm. We therefore expect that kaons are as (or even more) 
effective than electrons at screening the electric field in the CFL boundary layer. 



^If higher order corrections to the effective field theory for the Nambu-Goldstone bosons in the 
CFL phase were to make the 7r~ hghter than the K~ , we would of course obtain pion condensation 
rather than kaon condensation. 
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Figure 3: Particle number density profiles and electrostatic potential near a sharp 
interface between nuclear matter and CFL quark matter, including the possibility 
of kaon condensation on the CFL side. As in Fig. the transition occurs at 
fi = fin/^ = 365 MeV and an electron chemical potential /ig = 214 MeV. Solid 
curves show results which include effects due to strong interactions between baryons. 
Dashed curves show results where they are neglected. 



and should be included in a realistic description of the interface. This is easily 
done by including their contribution to the electric charge density in the source 
term of Poisson's equation. Particle profiles obtained upon including kaons in the 
CFL phase are shown Fig. ^. We see that the negative charge density in the CFL 
boundary layer is dominated by kaons, which are much more numerous than the 
electrons. This can be understood upon noting that wherever fif^ > rriK, the kaon 
density would rise without bound if the kaons were noninteracting bosons. It is 
only their Coulomb repulsion and the interactions encoded in the nonlinear effective 
Lagrangian for the kaons that stabilize their number density. Since electrons are 
fermions, their density is controlled by even in the absence of interactions. The 
figure also shows that the increased negative charge on the CFL side of the interface, 
due to the kaons, increases the electric field on the nuclear side of the interface, thus 
resulting in an even larger pileup of protons than we found in §^ Note that we 
have neglected higher derivative terms in the kaon effective field theory ([4.1|). This 



is justified wherever /Xg /A is small. The analysis of higher order terms in Ref. |T0 



suggests that they can be neglected for /if < (1 - 1.5) A. To see that this condition 



18 



-5 -4 -3 -2 -1 1 2 3 4 5 

e 

Figure 4: As in Fig. ^ but with number densities shown on a logarithmic scale 
and with the electrostatic potential shown as an inset. The solid curves are for 
iriK = 28 MeV, as obtained from ( [4.2| ) and as in Fig. |^. The dashed curves are for 
iriK = 100 MeV. The strong interactions among the nucleons are included. 



is satisfied for the profiles shown in Fig. note that because of the presence of the 
protons on the nuclear side of the interface, ^'^^ is at most ~ /ie/2 on the CFL side. 

We see from Fig. ^ that if the kaon mass is greater than that obtained from 
( |4.2| ), which is after all only quantitatively valid at asymptotically large densities, 
the maximum kaon density is somewhat reduced. More significantly, since kaon 
condensation occurs only where n'f > tjik, the kaon condensate does not extend 
as deep into the CFL phase for larger values of mx- Because we continue to take 
the electrons to be massless, we find that, as in Fig. 0, the electron density in the 
CFL phase has a power law tail, with kpe ~ I/2; at large z. For rriK = 28 MeV, the 
amplitude of this power law tail is greatly reduced, however, because most of the 
screening on the CFL side is now done by the kaon condensate. For larger values of 
rriK, there is less charge in the kaon condensate because its depth is reduced, and 
the amplitude of the power law tail of the electron density consequently increases, 
as the profiles become more similar to those of Fig. |[ 

In Fig. ^, we show the profiles at an interface in a theory in which the mixing 
angle which relates the ordinary photon to the Q-photon of the CFL phase is not 
small: we take coscto = 0.5. We see, for example, that the electron density is 
discontinuous because electrons have charge e on the nuclear side of the interface 
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e 

Figure 5: As in Fig. but for two values of the mixing angle uq. The solid curves are 
for costto = 1, as in all other figures. The dashed curves are for coscto = 0.5. The 
strong interactions between nucleons have been turned off. This figure illustrates 
the consequences of the fact that e < e and Eq > E hy dramatically exaggerating 
them. 

and charge e = ecosao on the CFL side. In nature ao ~ 1/20, so taking cosao = 1 
as we do in all other figures is certainly adequate. 

5 Surface Tension 

Having determined the profiles which characterize the minimal interface between 
CFL and nuclear matter, we now wish to investigate the stability of this sharp 
interface. We have already seen in Fig. |l] that if we neglect surface tension, the 
sharp interface must be less stable than a broad mixed phase region, denoted in 
that figure by the line AD. Before constructing this mixed phase, therefore, let us 
see what we can learn about the surface tension of the single sharp interface. 

There are two contributions to the surface tension. The first is that due to the 
Q CD- scale interface which we treat as sharp. We have no quantitative calculation of 
this surface tension, which we denote ctqcd, but we can estimate it by dimensional 
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analysis. The difference between the energy densities at 2; = 0^ and z = 0~ in, for 
example, the interface described by the solid curves in Fig. ^is about 325 MeV/fm^. 
If we estimate that the transition from CFL to nuclear matter occurs over a distance 
of about 1 fm, we expect 

c^QCD ~ 300 MeV/fm^ (5.1) 

This estimate is based only on dimensional analysis and could easily be a significant 
overestimate. 

The second contribution to the surface tension of the single interface is that due 
to the boundary layers whose profiles we have constructed. Again, we focus on the 
profiles shown as solid curves in Fig. ^. We must integrate the difference between 
the 2;-dependent energy density described by these number densities and the energy 
density of the undisturbed bulk CFL or nuclear phase. Solving the Poisson equation 
ensures that the pressure throughout the boundary layers is constant, and thus equal 
to that of the bulk phases. This means that the interface has no extra free energy. 
The energy density difference is therefore determined by the difference between the 
"/iA^" terms for the z-dependent profile and the undisturbed bulk phases. Hence, 
the contribution of the boundary layers to the surface tension is 

- r ( 

C"boundary layer / dz I /i^. [^n(^) ^n( ^)] 

J —00 V 

+ [fin - fJ'fiz)] np{z) - [fin - /Ue] ^p(-Oo) 
+ ljf{z)ne{z) - /Xene(-Oo) 



00 







+ / dz [ fif{z)n,iz)+i^f{z)nKiz) ] , (5.2) 



where ne^xiz) are the electron and kaon number densities in the CFL boundary layer 
and nn,p,e{z) are the neutron, proton and electron number densities in the nuclear 
boundary layer. We find 



CTboundary layer ^ 420 McV/fm (5.3) 

for the profiles in Fig. 0. Because /i„ and fip are much larger than /ig, and because 
Hn deviates little from its bulk value, the boundary layer contribution to the surface 
tension is in fact dominated by the contribution of the protons. We see that as a 
result of the development of charged boundary layers, and in particular as a result 
of the pileup of protons, the surface tension of the minimal CFL-nuclear interface 
is large and the dominant contribution is calculable. 

6 Mixed Phase 

If we neglect the surface tension and the Coulomb interaction then, as discussed in 
the introduction, the minimal CFL-nuclear interface is not stable because a broad 
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mixed phase region has lower free energy. The observation that a mixed phase may 
exist is neither new nor specific to the nuclear-CFL transition. In earher work, 
Glendenning showed that because of the existence of two independent chemical 
potentials, corresponding to conserved electric charge and baryon number, a mixed 
phase is a generic possibility wherever a first order phase transition occurs within 
a neutron star [Q. Subsequently, the first order transition between nuclear matter 



and unpaired quark matter was investigated by several authors [^, In 
these studies, wherein surface and Coulomb effects were either chosen to be small 
or neglected, a mixed phase was shown to be favored over a wide range of pressure. 

We begin this section with an explicit construction of the Glendenning-style 
mixed phase in the present context, neglecting the effects of surface tension and 
the Coulomb interaction. We shall restore these effects below. Phase coexistence 
is possible if the Gibbs condition of equal pressures and equal baryon chemical 
potentials is satisfied. (To maintain consistency of notation, we work with the 
quark chemical potential n; that for baryon number is 3fi.) At fixed fi, we find the 
electron chemical potential which satisfies the pressure equality condition 

-PcFL+kaons(/^, /^e) = -^nuclear (/^, /^e) , (6.1) 

where -PcFL+kaons describes the pressure of the negatively charged phase obtained by 
imposing /^e > on the CFL phase, thus creating a condensate of negatively charged 
kaons. This phase would have infinite (actually, nonextensive) free energy in bulk, 
due to the Coulomb repulsion. When it occurs in a mixed phase, the Coulomb energy 
cost (evaluated below) is finite. Note that no significant electric fields develop in 
the mixed phase, meaning that fi'f = /ig therein. As /Zg 7^ is uncompensated 
by any electrostatic potential, wherever fi^ > me we expect electrons and wherever 
fie > we expect a condensate of negatively charged kaons. 

The condition ( |6.1|) uniquely determines the pressure and function of 

the baryon chemical potential at which phase coexistence is possible. Clearly, the 
fie obtained as a solution to Eq. ( p.l|) does not correspond to the fie required to 
ensure the electrical neutrality of either phase. Electric charge neutrality is therefore 
enforced only as a global condition. This condition determines the volume fraction 
X of the CFL+kaon phase via 

XQcFL+kaons + (1 ^ X) Qnuclear = , (6.2) 

where Qnuciear IS the charge density of the nuclear phase, QcFL+kaons that of the CFL 
phase with kaons. 

In Fig. ^ we present the results of such a calculation. Several quantities that 
characterize: i) the neutral nuclear phase; ii) the nuclear-CFL mixed phase; and 
iii) the neutral CFL phase are shown. The electron chemical potential in the 
charge neutral nuclear phase is shown in solid black. It increases with increasing 
baryon chemical potential to compensate for electric charge density of the protons. 
The pressure of the electrically neutral nuclear matter phase is also shown. At 
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Figure 6: Behavior of the electron chemical potential and the pressure of 
homogeneous neutral nuclear and CFL matter and of the mixed phase, all as a 
function of the quark chemical potential fi. Only bulk free energy is included; 
surface and Coulomb energy is neglected. As in Fig. |l], the mixed phase occurs 
between A and D. The vertical line connecting B and C denotes the n at which 
the pressures of neutral CFL and nuclear matter are equal. This is where a sharp 
interface may occur. The pressure of the mixed phase exceeds that of neutral CFL 
or neutral nuclear matter between A and D. Were this the whole story, the mixed 
phase would evidently be favored over the sharp interface. 

H = 352 MeV, at the point A, the pressure of the electrically neutral nuclear 
matter phase and the negatively charged CFL+kaon phase coincide. Above this 
density, it becomes energetically favorable to construct a electrically neutral mixed 
phase wherein a positively charged nuclear phase coexists with a negatively charged 
CFL+kaon phase. At A, the negatively charged CFL phase containing kaons is 
constructed with the same value of fie as that in the charge neutral nuclear phase. 
For larger values of /i, the pressure and electron chemical potential of both the 
nuclear and CFL phases is obtained by solving Eq. ( |6.1| ). The pressure in the mixed 
phase is shown as the dot-dashed curve. In the mixed phase, fie decreases with 
increasing fi. The mixed phase ends when the electron chemical potential goes to 
zero, at the point D. At this point, fi = 408 MeV and the volume fraction of 
the CFL phase has reached x = 1- Beyond this density, we have homogeneous, 
electrically neutral, CFL matter. Fig. |^ also shows the critical point where the 
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pressure of electrically neutral nuclear matter equals that of Q-neutral CFL matter. 
This corresponds to the points labelled B and C, where fiBC = 365 MeV and 
fie = 214 MeV. This corresponds to the sharp interface, discussed in §|| and §H and 
shown in Fig. ^. We shall show below that this sharp interface is favored, even 
though it is evidently not favored by the bulk free energies, if it is stabilized by 
a sufficiently large surface tension. Note that in the neutral CFL matter to the 
right of the sharp interface, He is the same as that at B while /if = as a result 
of the electric field at the interface. In this way, the bulk CFL matter remains 
neutral (with no electrons or kaons) even in the presence of a large /ig. In the 
mixed phase, by contrast, there are no significant electric fields, /if = /te decreases 
continuously, and the CFL+kaon matter is negatively charged. The kaon condensate 
occurs only in that region of the mixed phase where /tg > ttlk- In Fig. |^, /tg = mx 
at /i = 398 MeV. Between this point and D, the negatively charged CFL component 
of the mixed phase contains electrons but no kaons. 

We must now evaluate the surface and Coulomb energy costs associated with 
the mixed phase, in order to make a fair comparison between it and the sharp 
interface. We have seen in §^ that the surface tension between a region of CFL 
matter and a region of nuclear matter is very large if the regions themselves extend 
over distances much greater than the electron Debye length Ae ~ 10 fm. From 
this we conclude that a mixed phase is only possible if the size of the CFL and 
nuclear regions within the mixed phase are comparable to or smaller than Ag. If 
this condition is satisfied, the particle number densities will be reasonably constant 



across a single CFL region or across a single nuclear region. (See Refs. [^, |3^, |38 
for demonstrations, in different contexts, that the relevant length scale is within a 
factor of two of Ag.) In a sense, if the characteristic length scale ro of the domains 
within the mixed phase is smaller than Ae, the mixed phase region is "all boundary 
layer" . If tq < Ae the surface tension whose free energetic cost we must yet take into 
account is dominated by o"qcd; that of the QCD-scale micro-boundaries between 
CFL and nuclear domains. If ro < Ae we may neglect corrections due to the (small) 
variation of particle densities within regions of size tq. 

For what value of cjqcd is the mixed phase energetically favored? To answer this, 
we will treat ctqcd as an independent parameter, constant across the whole mixed 
phase. In fact, ctqcd depends on both /t and /te,Qbut in our simple parametric study 
we will ignore this dependence and adopt the method outlined in Ref. to estimate 
the surface and Coulomb energy cost of the mixed phase. The mixed phase can 
be subdivided into electrically neutral unit cells called Wigner-Seitz cells, as in the 
analysis of the inner crust of a neutron star where droplets of charged nuclear matter 



coexist with a negatively charged fluid of neutrons and electrons [ 39 1 . In the present 



context, each Wigner-Seitz cell will contain some positively charged nuclear matter 



^The difference between tfie energy densities of the two pliases found at a given location within 
the mixed phase varies between about 300 and 350 MeV/fm^ as the mixed phase is traversed. 
As for the single sharp interface, therefore, naive dimensional analysis suggests a microboundary 
surface tension ctqcd ~ 300 MeV/fm^. 
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and some negatively charged CFL+kaon matter. Although at low temperature 
these unit cells form a Coulomb lattice, the interaction between adjacent cells can 
be neglected compared to the surface and Coulomb energy of each cell. In this 
Wigner-Seitz approximation, the surface and Coulomb energy per unit volume are 
fairly straightforward to calculate. They depend in general on the geometry and are 
given by 



E' = , (6.3) 

E^' = 27r aem/d(x) (AQf , (6.4) 

where d is the dimensionality of the structure {d = 1,2, and 3 correspond to Wigner- 
Seitz cells describing slab, rod and droplet configurations, respectively), a is the 
surface tension and AQ = Qnuciear — QcFL+kaons is the charge density contrast 
between the two phases. The other factors appearing in Eqs. ( |6.3|) , (|6.4|) are: x, 
the fraction of the rarer phase which is equal to x where x ^ 0.5 and (1 — x) where 
0.5 ^ X ^ 1; ro, the radius of the rarer phase (radius of drops or rods and half- 
thickness of slabs); and fd{x), the geometrical factor that arises in the calculation 
of the Coulomb energy which can be written as 

/.(x)-^ + . (6.5) 

The first step in the calculation is to evaluate tq by minimizing the sum of E'^ and 
E^ . The result is 

d X CTqcd 



We then use this value of tq in Eqs. ( |6.3| ) , ( |6.4| ) to evaluate the surface and Coulomb 
energy cost per unit volume 



E' + E^ = '- (47r d' Ux) x'Y" {AQf' . (6.7) 



We must now compare this cost to the bulk free energy benefit of the mixed phase. 

The lowest curve in Fig. ^ shows AQ, the difference between the free energy 
density of the mixed phase (calculated without including the surface and Coulomb 
energy cost) and the homogeneous electrically neutral nuclear and CFL phases 
separated by a single sharp interface. For fi ^ 365 MeV, Afl = Qmixcd — ^nuclear; for 
fi ^ 365 MeV, where ^cfl ^ ^nuclear, Afl = ilmixed — ^CFL- AQ is the difference 
between the Pnm/cfl and PMixed curves between A and D in Fig. The mixed 
phase has lower bulk free energy, so AQ, plotted in Fig. |^, is negative. 

The remaining curves in Fig. ^ show the sum of the bulk free energy difference 
AQ and {E^ + E'~^), the surface and Coulomb energy cost of the mixed phase 
calculated using Eq. (|6.7|) for droplets, rods and slabs for three different values of 
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Figure 7: The free energy difference between the mixed phase and the homogeneous 
neutral nuclear and CFL phases. In the lowest curve, the surface and Coulomb 
energy costs of the mixed phase are neglected, and the mixed phase therefore has 
the lower free energy. Other curves include surface and Coulomb energy for different 
values of ctqcd and different mixed phase geometry. As ctqcd increases, the surface 
and Coulomb price paid by the mixed phase increases. 



ctqcd- Careful inspection of the figure reveals that for any value of ctqcd, the mixed 
phase is described as a function of increasing density by a progression from drops 
to rods to slabs of CFL matter within nuclear matter to slabs to rods to drops of 
nuclear matter within CFL matter. This is the same progression of geometries seen 



in the inner crust of a neutron star or in the mixed phase at a first order phase 
transition between nuclear matter and a hadronic kaon condensate or unpaired 
quark matter [^. We have also checked that for ctqcd = 10 and 40 MeV/fm^, with 
the mixed phase geometry at any x taken to be that favored, the sizes of regions of 
both the rarer and more common phases (ro and its suitably defined counterpart) 
are always less than 5 — 6 fm. As this is less than Ae, we are justified in our neglect 
of any spatial variation within regions of a single phase, and are justified in our 
assertion that ctqcd is the relevant surface tension in the calculation of . 

For any given ctqcd? the mixed phase has lower free energy than homogeneous 
neutral CFL or nuclear matter wherever one of the curves in Fig. ^ for that ctqcd is 
negative. We see that much of the mixed phase will survive if ctqcd — 10 MeV/fm 
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while for ctqcd ^ 40 MeV/fm^ the mixed phase is not favored for any This 
means that if the QCD-scale surface tension (Tqcd ^ 40 MeV/fm^, the single sharp 
interface with its attendant boundary layers, described in previous sections, is free- 
energetically favored over the mixed phase. 



7 Looking Ahead to Neutron Star Structure and Collisions 

According to the considerations of this paper, neutron stars plausibly consist of 
nuclear matter of relatively low density floating on a dense CFL core, with a baryon 
and energy density discontinuity of about a factor of two at the astrophysically 
sharp interface separating them. Further work is certainly required to determine 
the location in /i, and ultimately in pressure and radius, of such an interface. A 
better understanding of the QCD-scale surface tension is required in order to confirm 
that the single sharp interface is stable against the formation of a broad mixed phase 
region. A mixed phase region would have distinctive characteristics. Its transport 
properties are very different from those of the uniform CFL state: neutrino mean free 
paths, which are long in the CFL phase |53], are very short in the mixed phase due to 



coherent scattering |^^. Ultimately, therefore, features in the temporal distribution 
of neutrinos emitted by a supernova and detected in an underground detector may 
allow a conclusive determination of whether a mixed phase region does or does not 
form. For the present, naive dimensional analysis suggests that ctqcd ~ 300 MeV, 
significantly greater than the minimum ctqcd ~ 40 MeV required to ensure that the 
sharp interface is favored. If present, the sharp interface has interesting properties 
at the tens of fermi length scale: a large pileup of protons on the nuclear side of the 
interface and a large amplitude kaon condensate on the CFL side of the interface. 

From the point of view of neutron star structure, though, the most important 
consequence of its presence is simply the discontinuous change in the density by 
about a factor of two. This will have quantitative effects on the mass vs. radius 
relationship for neutron stars with quark matter cores. It may also have qualitative 
effects on the gravitational wave profile emitted during the inspiral and merger 
of two compact stars of this type. One may expect characteristic features in the 
gravitational waves to arise both when the outer nuclear matter portions of the star 
begin to deform each other and then somewhat later when the denser CFL cores 
begin to deform. Neutron stars with two important length scales (the star radius 
and the core radius) will introduce features on two timescales in the gravitational 
waves produced late in an inspiral event. 
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